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The objective of this paper is to provide, for the problem of univariate symmetry (with respect 
to specified or unspecified location), a concept of optimality, and to construct tests achieving 
such optimality. This requires embedding symmetry into adequate families of asymmetric (local) 
alternatives. We construct such families by considering non-Gaussian generalizations of classical 
first-order Edgeworth expansions indexed by a measure of skewness such that (i) location, scale 
and skewness play well-separated roles (diagonality of the corresponding information matrices) 
and (ii) the classical tests based on the Pearson-Fisher coefficient of skewness are optimal in 
the vicinity of Gaussian densities. 
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1. Introduction 

1.1. Testing for symmetry 

Symmetry is one of the most important and fundamental structural assumptions in 
statistics, playing a major role, for instance, in the identifiability of location or intercept 
under nonparamctric conditions: see [2, 14, 15]. This importance explains the huge variety 
of existing testing procedures of the null hypothesis of symmetry in an i.i.d. sample 
Xi,..., X n ; see [8] for a survey. 

Traditional tests of the null hypothesis of symmetry - the hypothesis under which X\ — 

9 = —{X\ — 9) for some location 6 € M, with = standing for equality in distribution - are 

based on standardized empirical third-order moments. Let m^ l \9) := nT 1 Y^i=\{Xi — 9) k 

and := m^\x^), where X'"' := nT 1 J27=i ^i- When the location 9 is specified, 
the test statistic is 

S^\e):=n^mt\e)/(m^\9)f 2 , (1.1) 
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the null distribution of which, under finite sixth-order moments, is asymptotically stan- 
dard normal. When 9 is unspecified, the classical test is based on the empirical coefficient 
of skewness 

b^:=m^/sl (1.2) 

where s n := (mj"') 1 '' 2 stands for the empirical standard error in a sample of size n. More 
precisely, this test relies on the asymptotic standard normal distribution (still under finite 
moments of order six) of 

Si n) := n^m^/rf 5 " 6s n m ^ + 9a n) 1/2 , (1-3) 

which, under Gaussian densities, asymptotically reduces to \/n/Qo{ . 

These two tests are generally considered as Gaussian procedures, although they do 
not require any Gaussian assumptions and despite the fact that none of them can be 
considered optimal in any Gaussian sense, since asymmetric alternatives clearly cannot 
belong to a Gaussian universe. Despite the long history of the problem, the optimality 
features of those classical procedures thus are all but clear, and optimality issues, in that 
fundamental problem, remain essentially unexplored. 

The main objective of this paper is to provide this classical testing problem with 
a concept of optimality that confirms practitioners' intuition (i.e., justifying the b^- 
based Gaussian practice) , and to construct tests achieving such optimality. This requires 
embedding the null hypothesis of symmetry into adequate families of asymmetric alter- 
natives. We therefore define local (in the LeCam sense) alternatives indexed by location, 
scale and a measure of skewness in such a way that: 

(i) Location, scale, and skewness play well-separated roles (diagonality of the corre- 
sponding information matrices). 

(ii) The traditional tests based on b± (more precisely, based on S% given in (1.3)) 
become locally and asymptotically optimal in the vicinity of Gaussian densities. 

As we shall see, part (ii) of this objective is achieved by considering local first-order 
Edgeworth approximations of the form 

<j>(x -9) + n- 1/2 £{x - 6)<j>(x - 9){{x - 9) 2 - «), (1.4) 

where <fi as usual stands for the standard normal density, k (= 3) is the Gaussian kurtosis 
coefficient, 9 is a location parameter, and £ a measure of skewness. Adequate modifica- 
tions of (1.4), playing similar roles in the vicinity of non-Gaussian standardized symmetric 
reference densities /i, are proposed in (2.2). 

The resulting tests of symmetry (for specified as well as for unspecified location 9) 
are valid under a broad class of symmetric densities, and paramctrically efficient at 
the reference (standardized) density f\. Of particular interest are the pseudo- Gaussian 
tests (associated with a Gaussian reference density; see Proposition 3.6), which require 
finite moments of order six and appear to be asymptotically equivalent (under their 
specified-^ version as well as under the unspecified- 9 one) to the test (1.3) based on b^"\ 
and the Laplace tests (associated with a double-exponential reference density; see Propo- 
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sition 3.6), which only require moments of order four and are closely related with the 
tests against Fechner asymmetry derived in [3]. 

These tests are of a parametric nature. Since the null hypothesis of symmetry enjoys 
a rich group invariance structure, classical maximal invariance arguments naturally bring 
signs and signed ranks into the picture. Such a nonparametric approach is adopted in 
a companion paper [4], where we construct signed-rank versions of the parametrically 
efficient tests proposed here. These signed-rank tests are distribution-free (asymptotically 
so in case of an unspecified location 9) under the null hypothesis of symmetry, and 
therefore remain valid under much milder distributional assumptions (for the specified 
location case, they are valid in the absence of any distributional assumption). 

The main technical tool throughout the paper is LeCam's asymptotic theory of statisti- 
cal experiments and the properties of locally asymptotically normal (LAN) families. LAN 
has become a standard tool in asymptotics: sec [12] or Chapters 6-9 of [17] for details. 
Log-likclihoods in a LAN family with fc-dimcnsional parameter # admit local quadratic 
approximations of the form t'A^ 1 ' — ^t'T^t, where the random vector A^ 1 ', called 
a central sequence, is asymptotically normal Af(T^T,T^) under sequences of parameter 
values of the form i9 + n~ x / 2 T [local alternatives). Let </>*(A) be an optimal test (uni- 
formly most powerful, maximin, most stringent, . . . ) in the Gaussian shift model describ- 
ing a hypothetical observation A with distribution in the family {J\f(T#T, r#) | r € R fe } 
(r# specified). Those families are extremely simple, and optimal tests in that context are 
well known (see, e.g., Section 11.9 of [11]). Then, the sequence (f)*(A^) is a sequence of 
locally asymptotically optimal (locally asymptotically uniformly most powerful, maximin, 
most stringent, . . . ) tests for the original problem - where optimality is based on the local 
convergence of risk functions to the risk functions of Gaussian shift experiments. More 
analytical characterizations can be found, for instance, in [5]. 

1.2. Outline of the paper 

The problem we are considering throughout is that of testing the null hypothesis of 
symmetry In the notation of Section 1.1, £ (see (2.2) for a more precise definition) is thus 
the parameter of interest; the location 9 and the standardized null symmetric density f\ 
either are specified or play the role of nuisance parameters, whereas the scale a (not 
necessarily a standard error) always is a nuisance. 

The paper is organized as follows. In Section 2.1 we describe the Edgeworth-type fam- 
ilies of local alternatives, extending (1.4), that we are considering. Section 2.2 establishes 
the local and asymptotic normality (with respect to location, scale and the asymme- 
try parameters) result to be used throughout; actually, we establish a slightly stronger 
version of LAN, called ULAN (uniform LAN), which allows us to handle the problems 
related with estimated nuisance parameters. The classical LeCam theory then is used 
in Section 3.1 for developing asymptotically optimal procedures for testing symmetry 
(£ = 0), with specified or unspecified location 9 but specified standardized symmetric 
density fx. The more realistic case of an unspecified f\ is treated in Section 3.2, where 
we obtain versions of the optimal (at given /i) tests that remain valid under gi ^ f\, 
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for specified (Section 3.2.1) and unspecified (Section 3.2.2) location 6, respectively. The 
particular case of pseudo-Gaussian procedures (which are optimal for Gaussian /i but 
valid under any symmetric density with finite moments of order six) is studied in detail in 
Section 3.3 and their relation with classical tests of symmetry is discussed. We also show 
that the Laplace tests (which are optimal for double-exponential /i but valid under any 
symmetric density with finite fourth-order moment) are closely related to the Fechncr- 
type tests derived in [3] . The finite-sample performances of these tests are investigated 
via simulations in Section 4, where they are applied to the classical skew-normal and 
skew-t densities. 



2. A class of locally asymptotically normal families of 
asymmetric distributions 

2.1. Families of asymmetric densities based on Edgeworth 
approximations 

Denote by X<"> := (x[ n) , . . . , A^ n) ), n G N an i.i.d. n-tuple of observations with common 
density /. The null hypotheses we are interested in are: 

(a) The hypothesis Hg 1 ^ of symmetry with respect to specified location 6 <G R: un- 

(n) 

der H g , the Xi's have density function 

x^f(x):=a- 1 h((x-6)/a) (2.1) 

(all densities are over the real line, with respect to the Lebesgue measure) for 
some unspecified a € Rj , where fi belongs to the class of standardized symmetric 
densities 

To := : h\{— z) = h\(z) for all z G M and J h\(z) dz = 0.75 

The scale parameter a (associated with the symmetric density /) that we are 
considering here thus is not the standard deviation, but the median of the absolute 
deviations \Xi — 9\; this avoids making any moment assumptions. 

(b) The hypothesis := UeeR^e"^ °f symmetry with respect to unspecified loca- 
tion (and scale): there exist (8, a) such that the X^s have density (2.1). 

In both cases, the standardized density fi may be specified (Section 3.1) or not (Sec- 
tions 3.2-3.4). The specified- f± problem, however, mainly serves as a preparation for the 
more realistic unspecified- /i one. 

As explained in the introduction, a characterization of efficient testing requires the 
definition of families of asymmetric alternatives exhibiting some adequate structure, such 
as local asymptotic normality, at the null. For a selected class of densities / enjoying 
the required regularity assumptions, we therefore are embedding the null hypothesis of 
symmetry into families of distributions indexed by 9 £ R (location), a G Rj (scale) and 
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a parameter £ £ M characterizing asymmetry. More precisely, consider the class T\ of 
densities f\ satisfying: 

(i) (symmetry and standardization) /i € Jo; 

(ii) (absolute continuity) there exists f\ such that, for all z\ < Zi, 

h(z2)-fi(zi)= r fi(z)dz; 

(iii) (strong unimodality) z i— > 0^ (z) := —j\(£)Jf\(z) is monotone increasing; 

(iv) (finite Fisher information) /C(/i) := z 4 ^ (z)/i(z) dz, hence also, under 
strong unimodality, 

/+oo /-+oo 
2 /i (z)/ 1 (z)dz and J(fi):= z 2 ^ (z)A(z) dz 

-oo J —oo 

are finite; 

(v) (polynomial tails) 

/•OO 

/ f 1 {z)dz = 0(y- fi ) asy^oo 
•'s/ 

for some /3 > and 

<?!>/! (z) =o(z /3/2 " 2 ) asz^oo. 

That class J-\ thus consists of all symmetric standardized densities f\ that are absolutely 
continuous, strongly unimodal (that is, log-concave) and have finite information X(/i) 
and v7(/i) for location and scale, and, as we shall see, IC(fi) for asymmetry, with tails 
satisfying (v). 

For all /i £ Ji, denote by k(/i) := J(fi)/I(f{) the ratio of information for scale and 
information for location; k(/i), as we shall sec, for Gaussian density (/i = 4>\) reduces to 
kurtosis (k(c/>i) = 3), and can be interpreted as a generalized kurtosis coefficient. Finally, 
write Pg 1 ^ ^ for the probability distribution of X'™' when the Xi's are i.i.d. with density 

/(*) = °- x h (^) - ^- X h (^) ( (^) 2 - «(/i)) /D* - *l < 

- sig^Oa" 1 /! f— J {/[a - > sign(-£)a|z*|] - /[x - < S ign(£>|**|]}. 



Here 9 € K and cr g R + clearly are location and scale parameters, £ € K is a measure 
of skewness, (strictly positive for f\ e T\) the generalized kurtosis coefficient just 

defined and z* the unique (for £ small enough; unicity follows from the monotonicity 
of 4> h ) solution of /i(z*) = £/i(z*)((z*) 2 - k(/i))- The function / defined in (2.2) is 
indeed a probability density (non- negative, integrating up to one), since it is obtained 
by adding and subtracting the same probability mass 
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on both sides of 9 (according to the sign of £). Note that £ > implies f(x) = for 
x - 9 < -a\z*\ and f(x) = 2a~ 1 h((x - 6) /a) for x — 6 > cr\z*\. Moreover, x i-» f(x) is 
continuous whenever fi(x) is; vanishes for x < 9 + az* if £ > 0, for x > 9 + az* if £ < 0; 
and is left- or right-skewed according as £ < or £ > 0. As for z* , it tends to — oo as £ J, 0, 
to oo as £ 1 0; in the Gaussian case, it is easy to check that \z*\ = 0(|£| -1 / 3 ) as £ — > 0. 

The intuition behind this class of alternatives is that, in the Gaussian case, (2.2), with 
£ = n _1 ' 2 T yields (for x <G [9±az*]) the first-order Edgeworth development of the density 
of the standardized mean of an i.i.d. n-tuple of variables with third-order moment 6tct 3 
(where standardization is based on the median a of absolute deviations from 9). For 
a "small" value of the asymmetry parameter £, of the form n -1 / 2 r, (2.2) thus describes 
the type of deviation from symmetry that corresponds to the classical central limit con- 
text. Hence, if a Gaussian density is justified as resulting from the additive combination 
of a large number of small independent symmetric shocks, the locally asymmetric / re- 
sults from the same additive combination of independent but slightly skew shocks. As we 
shall see, the locally optimal test in such a case happens to be the traditional test based 
on (see (1.2)). 

Besides the Gaussian one (with standardized density (j>i(z) := y / a/27texp(— az 2 /2)), 
interesting special cases of (2.2) are obtained in the vicinity of: 

(i) The double-exponential or Laplace distributions, with standardized density 

h(z) = f c (z):=(l/2d)eM-\z\/d), 

1(h) = l/d 2 , J{h) = 2 and JC(h) = 24d 2 . 

(ii) The logistic distributions, with standardized density 

fi(z) = ho g (z) := Vbexp(-Vbz)/(l + exp(-^z)) 2 , 

Ah) = b/3, J{fi) = (12 + 7t 2 )/9 and K{h) = 7t 2 (120 + 7n 2 )/A5b. 

(iii) The power-exponential distributions, with standardized densities 

fi(z) = fex Pv (z) := C cxPv exp(-(^z) 2 "), 

V G N , 2(h) = 2g 2 V T(2 - 1/2 V )/T(1 + 1/2 V ), J(h) = l + 2 V and JC(h) = 2g v r,/ 
T(l + l/27j) (t ne positive constants C eX p^, a, b, d and g v are such that f\ e ^i). 

Although not strongly unimodal, the Student distributions with v > 2 degrees of free- 
dom also can be considered here (strong unimodality indeed is essentially used as a suffi- 
cient condition for the existence of z* in (2.2) - an existence that can be checked directly 
in the Student case). Standardized Student densities take the form 

h(z) = f tv (z) := C tv (1 + a^/z/r^ 1 )/ 2 , 

with 1(h) = a v (v + l)/(y + 3), J(h) = 3(f + l)/(v + 3) and K.(h) = I5v(v + l)/a v (v- 
2)(z^ + 3) (Ct v and a v are normalizing constants). Note that the corresponding Gaussian 
values, namely l(4>i) = a ps 0.4549, J(4>i) = 3 and JC(4>i) = 15/a, are obtained by taking 
limits as v — > 00. 
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Figures 1 and 2 provide graphical representations of some densities in the Gaussian 
(fi = 4>i) and double-exponential (/i = fc) Edgcworth families (2.2), respectively. In the 
Gaussian case, the skewed densities are continuous, while the double-exponential ones, 
due to the discontinuity of fc(%) at x = 0, exhibit a discontinuity at the origin. 



2.2. Uniform local asymptotic normality (ULAN) 

The main technical tool in our derivation of optimal tests is the uniform local asymptotic 
normality (ULAN), with respect to # := (0,a,£)', at (6, a, 0)', of the parametric families 

v h ■= U v i% ■= U i p o n Uh I ^ R,C e M}, (2.3) 

(T>0 (T>0 

where fi E T\. More precisely, the following result holds (see the Appendix for proof). 

Proposition 2.1 (ULAN). For any fx € T\, 8 G K, and a £ the family is 

ULAN at {6,a,0y, with (writing Z % for z\ n \e,a) := cr^P^"' - 6) and <f> fl for -fi/fl) 
central sequence 





Figure 2. Graphical representation of the double-exponential Edgeworth family (2.2) (/i = fc), 
for £ = 0, 0.05, 0.10 and 0.15. 
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and full-rank information matrix 

- 2 Ah) 
o 

where 7 (/i) := AC(/i) - J 2 {h)/I{h). 

More precisely, for any := {8^ , crW , 0)' such that 0<» - (9 = 0(n^ 1 / 2 ) and crW - 
a = 0(n~ 1 / 2 ) , and for any bounded sequence t^ 1 ' = (f( n ) , s^"\t^)' gR 3 , we Ziave, under 

dP (n) 



»(n) , / i?("'+» 1 - 1 / 2 t(");/i 

A ^(») +n -l/2 x („) /l? (n) ;/l — 10 Sl („) 



X («)/ A («) (l9 («) ) _ 1 (n)' r (^) T W +Qp (l) 



and 



The diagonal form of the information matrix T ^ (#) confirms that location, scale and 
skewness, in the parametric family (2.3), play distinct and well-separated roles. The 
practical consequences of that orthogonality (in the sense of information) property, as 
we shall see, are twofold: 

(a) The fact that location and scale are unspecified has no cost in terms of efficiency 
and under specified standardized density f\ when testing for symmetry. 

(b) Substituting root-n-consistent (and, in principle, duly discretized: see assumption 
(C2) below) estimators for the true values has no impact on the asymptotic validity 
and local powers, under specified standardized density f\, of tests for symmetry. 

Note that orthogonality between the scale and skewness components of A.^ ($) automat- 
ically follows from the symmetry of /i, while for location and skewness, this orthogonality 
is a consequence of the definition of The Gaussian versions of (2.4) and (2.5) are 



and 

respectively (recall that a « 0.4549). 





Optimal tests for symmetry 



1071 



3. Locally asymptotically optimal tests 

The various test statistics described in this section are listed, for easy reference, in Table 3, 
placed at the end of Section 5. 

3.1. Locally asymptotically optimal tests: Specified density 

For specified /i S T\ , consider the null hypothesis 

: = U { P ^o ;/l } 

of symmetry with respect to some specified location 9, and the null hypothesis 

K ] = U U { P £o;/J 

of symmetry with respect to unspecified 9. ULAN and the diagonal structure of (2.5) im- 
ply that substituting discretized root-rt-consistcnt estimators 9 and a for the unknown 9 
and a has no influence, asymptotically, on the £-part of the central sequence. 

Recall that a sequence of estimators A*-™) defined in a sequence of experiments {P^ | 
A € A} indexed by some parameter A is root-n- consistent and asymptotically discrete if, 
under P^ 1 , as n — > oo: 

(CI) A(")-A = P (n- 1 /2). 

(C2) The number of possible values of X^ in balls with ©(n" 1 / 2 ) radius centered at A 
is bounded as n — > oo . 

An estimator X^ satisfying (CI) but not (C2) is easily discretized by letting, for some 
arbitrary constant c > 0, A^ := (cn 1 / 2 )" 1 sign(A'"- ) )[cn 1 / 2 |A^"- ) |] , which satisfies both 

(CI) and (C2). Subscripts # in the sequel are used for estimators (6#, <5#, ■ ■■) satisfying 
(CI) and (C2). It should be noted, however, that (C2) has no implications in practice, 
where n is fixed, as the discretization constant c can be chosen arbitrarily large. 

It follows from the diagonal form of the information matrix (2.5) that locally uni- 
formly asymptotically most powerful tests of (resp., of ) can be based on 

Ag 3 (0,<r # ,O) (resp., on A^ 3 (<9 # , <r # , 0)), hence on T { f f(9,a # ) (resp., on (§#,&#)), 
where 

n 

t£\0,<t) := = J2<t>fAZi((>,<r))(Zi(0,cr) - K (h))- (3-1) 

Root-n-consistent (under the null hypothesis of symmetry) estimators of 9 and a that do 
not require any moment assumptions are, for instance, the medians 9 := Med(x|™') and 
a := Med(|X} n ^ — 6\) of the X^'s and of their absolute deviations from 9, respectively. 
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The following proposition then results from classical results on ULAN families (see, 
e.g., Chapter 11 of [11]). 

Proposition 3.1. Let /i £ Ji. Then: 

(i) Tf™\§#,a#) =Tj™'(6,o~) +op(l) is asymptotically normal, with mean zero under 

Pg^o-fi' mean T 7 1 ^ 2 (/i) un d er ^i} l ] rn -i/2 T .f 1 an d variance one under both. 

(ii) The sequence of tests rejecting the null hypothesis of symmetry (with standard- 
ized density f\) whenever Tj™ {6,a#) (resp., Tj™ (0#, exceeds the (1 — a) 
standard normal quantile z a is locally asymptotically most powerful at asymp- 
totic level a for (resp., for H^' ) against (J ?>0 U CTe R+{ P e7<7,€;/i} ( res P-, 

Uf >o UeeK Uo-eR+ $,<?,£ ;/i })■ 

This confirms that unspecified location 9 and scale a do not induce any loss of efficiency 
when the standardized density f\ itself is specified. 
The Gaussian version of (3.1) is 




thanks to the linearity of Gaussian scores, it easily follows from a traditional Slutsky 
argument that 9 and a in T^(9,a) need not be discretized. Under Gaussian den- 
sities, both T^\§,a) and T^'(6,a) are asymptotically equivalent to Tfj™'(X( n \ a) = 

(na 3 /6) 1/2 TO3 Il) /o- 3 = y^n~/6 b[ n) +o P (l), that is, to S^ n) given in (1.3). The latter is thus 
locally asymptotically optimal under Gaussian assumptions, whether 9 is specified or not, 
whereas the specified-^ test based on m < 3 l) {9)/(m i Q l) {9)) 1 / 2 (more precisely, on s[ n) {9) 
given in (1.1)) is suboptimal. The fact that rrij" (0) yields a better performance than 
(9) under specified location 9 (see the comments after Proposition 3.5 for a com- 
parison of local powers) looks puzzling at first sight. The reason is that orthogonality, in 
the Fisher information sense, between asymmetry and location, is a "built-in" feature of 
Edgeworth families. Contrary to S^, which is shift invariant, m^ L \9) and s[ n \9) are 

sensitive to location shifts. Therefore, tests based on s[ (6) are "wasting away" some 
power on location alternatives (which are irrelevant when 9 is specified), to the detriment 
of asymmetry alternatives. 

Locally asymptotically maximin two-sided tests are easily derived along the same lines. 

3.2. Locally asymptotically optimal tests: Unspecified density 

The tests based on (3.1) achieve local and asymptotic optimality at correctly specified /i, 
which sets the parametric efficiency bounds for the problem, but has limited practical 
value, as these tests are not valid anymore under density g\ ^ f\. If Proposition 3.1 is to 
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be adapted to the more realistic null hypotheses := |J gi and H {n) := \J gi ?4? } 

under which the (symmetric) density remains unspecified, the test statistic Tj™ needs to 
be adapted in order to cope with three problems: its centering under the null and g\ 7^ f\ 
iTf™' in the previous section only had to be centered under density g\ = its scaling 
(same remark) and (still under the null and 31 7^ /1) the impact of the substitution of 
estimators a (and 9) for the unspecified values of a (and 8) on its asymptotic distribution. 

3.2.1. Specified location 

in) 

Let us first assume that both 9 and a are specified. The test statistic in Propo- 
sition 3.1 is essentially a scaled version of Aj™. 3 . More generally, write Ai™. 3 (/«) for 
nT 1 / 2 X]ILi 0/i ~ K )i where K <E denotes a strictly positive real number; an ad- 

equate sample-based value will be selected later on. Note that A^ 3 (k) remains centered 

under Pg™^ . gi , irrespective of the choice of k. Indeed, the functions z 1— > <f)fi ( z ) z2 an d z <— > 
4>f 1 (z) are skew symmetric, and their expectations under any symmetric density arc auto- 
matically zero - provided that they exist. The variance under Pgl o- gi °f ^/™3( K ) 1S then 



l gi (h) ■= E 31 U h (Z t )(Z? - n)f] = fC gi (A) - 2nJ gi (A) + n 2 l gi (ft), 



where 



/DO rOO 
4,%(z) gi (z)dz, JM):= z 2 <Pl(z) gi (z)dz 
-00 J — 00 

and (still, provided that those integrals exist) 

r+00 

^3i(/i) : =/ z i 4> 2 fl {z)g 1 (z)Az. 



We know from LeCam's third lemma that, under Pglo- gi , the impact on A^' 3 (k) 
of an estimated scale depends on the asymptotic joint distribution (under Pg™^o 9l ) °f 

A^ 3 (k) and A^ 2 . More precisely, LeCam's third lemma (see, e.g., page 90 of [17]) tells 

(n) 

0,a,O; gi ' 



us that if, under P*-™ 



/ aS&(«) \ 

1°S — 7~\ ; 




as n->oo, (3.2) 



then, under Pg™^ +n _i /2T , A^ 3 (k) — > Af(fii + dn,df). In the present context, Propo- 
sition 2.1 yields that 

log 9 ;V 2 "° ;91 = tA ^ («. 0) - \r 2 a- 2 (J( 91 ) 1) + o P (l) (3.3) 
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as n — > oo, under P^ Q ; hence (still as n — > oo under P^ o-gi)' 



^(( J -o (A, )'(*. 4 



(3.4) 



with d| = cr ~ 2 («^(.9i) — 1) an d dn— to~ x (j)^ (z)(z 2 — n)(4> gi (z)z — l)gi(z) dz = (as 
the integral of a skew-symmetric function). We conclude that (3.2) holds with du = 0. 

LeCam's third lemma therefore shows that the effect on the asymptotic distribution of 
Aj™? 3 (k) of a 0(r7, -1 / 2 ) perturbation of a is asymptotically nil; the asymptotic linearity 
result of Proposition A.l and a classical argument on asymptotically discrete estima- 
tors (see, e.g., Lemma 4.4 in [10]) allow for extending this conclusion to the stochastic 
Op(n -1 / 2 ) perturbations induced by substituting a duly discrctized root-n-consistent 

~ (n) 

estimator a\ for a . Such a substitution consequently does not affect the asymptotic 
behavior of Ay™. 3 (re). 

For fi G T\ and gi G Tfa := {.<?i G F\:K, gi {fi) < oo} (due to strong unimodality, 
fcgi (/i) < oo also implies T gi (fi) < oo and J gi (/i) < oo), let 

7 (r °(/i)=7 (n) (/i,0,O :=^(/ 1 )-2 K (/ 1 )jW(/ 1 )+« 2 (/ 1 )jW(/ 1 ), (3.5) 

where 

n 

= X (n) flj a) := n -l £ ^ (z . (ffj (3 g) 

n 

J(n) (/l ) = J(n) (/i ^ ct) ;= n -l £ Z 2 {z . ^ a)) (3 ?) 



and 



£(») = £(») (/i ^ ^ a) ;= n -x J- Zi(6, a)cj>l (Zi(0, a)) (3.8) 



under Pg™^ arc consistent estimates of X gi (fi), J gi (/i) and JC gi (fx), respectively. Now 
in practice, 

Z {n) (fi), J [n) {h) and hence 7 (n) (/i)> cannot be computed from 

the observations and Zj(0, <7#) is to be substituted for Zi(8,a) in (3.6)-(3.8), yielding 
7^ n H/i A<5#)- This substitution in general requires a slight reinforcement of regular- 
ity assumptions. Along the same lines as above (LeCam's third lemma and asymptotic 
linearity), we easily obtain that 7 ( -™' ) (/i,#,o'#) — ^ n \fi, 0, cr) is op(l) under P^o-gi 

provided that the asymptotic covariance of 7<")(/i,6»,ct) and A g % is finite. A simple 
computation (and the strong unimodality of /i and g\) shows that a sufficient condition 
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for this is 



gx € F* fl := j/ii e ^ : y z 5 ^ l (z)< ? !), I1 (z)/ii(z)dz<oo|. 



(3.9) 



Defining the test statistic 



and the cross-information quantities 

/+oo 
0/i( z )09i( z )5i(^)dz, 
-oo 

/+oo 
z 2 4> h {z)(t) gi {z)g 1 {z)dz 
-oo 

and 

/+oo 
z 4 /l (z)0 gi (z) 5 i(z)dz 
-oo 

(which for fx £ J 7 ! and 91 € J 7 ^ are finite because of Cauchy-Schwarz), we thus have the 
following result. 

Lemma 3.1. Let f\ E Fx and gx E J 7 ^ ■ Then: 

(i) Tf™\0,a#) = Tf™'(9,a) +op(l) is asymptotically normal, with mean zero under 
F e,v,om> mean 



T 



£gi(fl,9l) ~ Jgi(fl,9l)Wl) +«(ffi)) +I gi (fx,gi)K(fx)K(gx) 



[JCgAh) - XTMMfi) +X fll (/i)« 2 (/i)] 1 / 



2 



(3.11) 



under , and variance one under both. 

(ii) The sequence of tests rejecting the null hypothesis '■= U f r» °f sym- 

metry with respect to specified 9 whenever T^\6, a #) exceeds the (1 — a) standard 
normal quantile z a is locally uniformly asymptotically most powerful at asymptotic 
level a for H { g l) against (J ?>0 U a eU+ { P e?2,€;/ J • 

The tests based on T{™'(6,a#) enjoy all the validity (under %g ') and optimality 

(against Uf>o UcreR+iPe"! £-/i}) properties one can expect. However, a closer look reveals 
that they are quite unsatisfactory on one count: under gx 7^ fx, their behavior strongly 
depends on the arbitrary choice of the concept of scale (here, the median of absolute 
deviations) . 
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Consider, for example, the Gaussian version of (3.10), which takes the form 

n 

where 

7 (")(^) = 7 (™)(0 1 ,0,a) =a 2 a- 6 774 n) (0) -6aa- 4 mi n) (0) + 9(T- 2 774" ) (0). 

The test based on T^\o,&) (here again, Slutsky's lemma allows for not discretizing <x) 
is a pseudo-Gaussian test, hence optimal under Gaussian assumptions; the asymptotic 
shift (3.11) is r^jQ/a under ~Pg^ n -i/i T .j, : and 

r[5a^ 4 (gi) - (9 + 3aK(gi))^ 2 (gi) +3/«(gi)][a 2 /i 6 (5i) - 6a/j 4 (.gi) + 9^2(51 )]~ 1/2 , 

where fJ-k(gi) '■= zk 9i( z ) dz, under n -i/2 T . gi ■ This asymptotic shift strongly de- 
pends on a, hence on our (arbitrary) choice of a scale parameter. Setting to one the 
standard deviation instead of the median of absolute deviations would significantly mod- 
ify the local behaviour of T^™' (9, a#) as soon as g\ 7^/1. This does not affect optimality 
properties (which hold under /1), but is highly undesirable. 

Now, that unpleasant feature of is entirely due to the choice of k = ft(/i) as 
a (nonrandom) centering in (3.10). That choice was entirely motivated by asymptotic 

(n) 

orthogonality considerations under P g ^o-/i> and does not affect the validity of the test. 
It follows that replacing k(/i) with any data-dependent sequence kS 71 ' such that KS n > — 
K (fi) =op(l) under Pg^o-A asymptotically has no impact on T^\o,a) under Pg™ ^ a .f x - 
Let us show that this sequence «W can be chosen in order to cancel the unpleasant 
dependence of the test statistic on the definition of scale. 
Provided that 

fi € J~i := {hi e Ji : z h> 4>h 1 (z) is differentiable, with derivative 
integration by parts yields 



J 9i(/i>5i) = / <t>h(z)9i(z)&z 

J — OO 

and 

Jgi(fi,9i) = 2 z4> fl (z)gi(z)dz+ / z 2 4> fl (z)g 1 (z)dz. 



Therefore, Z 9l (/i, 51), J gi (fi,9i) and re ai (/i,5i) := J~ Sl (/1, (/i,.9i) under Pg™i i0;gi 
are consistently estimated by 

n 

I W °(/0 =I ( " )o (/i,0,a) := - V^Z^a)), 

i=l 

2 " 1 ™ 

J- (n)o (/i) = J- (n)o (/i, 9, a) := - V Zi{9, a)<t» h (Z^O, a)) + - V Zf(9, a)<j> fl (Z t (9,<j)) 



n 

i=l 
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and 

K W °(/i) = <s Wo (/i,8 I ff):= 1 7 Wo (/i)/I (n)8 (/i), (3.12) 

respectively. 

Clearly, K^ n '°(fi) satisfies the requirement that K^ n '°(fi_) — k(/i) = op(l) under 
Pg^ o-.fr ^ n P rac tice, however, K^°(fi,6,a) cannot be computed from the observations, 
and n( n > o (fi,0,(T#), where Zi(0,a#) has been substituted for Zi(6,a), is to be used in- 
stead. As in the evaluation of 7^(/i) above (see (3.5)), this substitution requires mild 
additional regularity conditions. LeCam's third lemma then applies exactly along the 
same lines, implying that 

K^ {f u 6,a # )-K^ {h,e,a) = o v {l) 

under Pg^o-gi as soon as the asymptotic covarianccs of and J^°{f\) with 

A^™? 2 are finite. A simple computation (and the strong unimodality of fi and g\) shows 
that a sufficient conditions for this is 

,9i G T° h := G T* fl : j z 3 4) fl (z)4) hl (z)h 1 (z)dz < oo 

(3.13) 

and / ztj)f 1 (z)(f)/ ll (z)hi(z)dz < < 

J —oo 

(no redundancy, since 4>f t is not necessarily monotone). 

Emphasize the dependence of A^ 3 (k) on 8 and a by writing A^". 3 (re, 9, a): it follows 
from Lemma A. 5 in the Appendix that, for f\ 6 J-° and g\ G T ^ , the difference between 

AS3,(/«< n > o (/i,0,*#),0,fr#) and A^3( Kffl (/ ljffl ),^a) is o P (l) under P^o- Sl - Lettin § 
(still for fiG^i) 

T^°(6,a):= ^ V0 A (^(M)(^(^) - ^(/t)), (3.14) 

where 

7 (n) °(/i) - 7 (n) °(/i,e,^) := K (n) (/i) - 2 K ( " )o (/ 1 )^ ( " ) (/ 1 ) + (« (n)o (/i)) 2 I (n) (/i), 
we thus have the following result. 

Proposition 3.2. Le£ /i £ J-^ and pi 6 J 7 ^ . Then: 



(i) Tf™*°(9,a#) =Tj™'°(6, cr) + op(l) zs asymptotically normal, with mean zero under 

£gi(h,9l)-Jgdfl>9l)Kgi(h,9l) 



T>( n ) 



[/C S1 (/i) - 2J S1 (A, 5 i) (A, 



(3.15) 



under Py' -mi an d variance one under both. 
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(ii) The sequence of tests rejecting the null hypothesis Tig := (J _ „ T~iq „ of sym- 
metry (with specified location 9, unspecified scale a and unspecified standardized 
density g i 6 J 7 ]^) whenever (6, a#) exceeds the (1 — a) standard normal quan- 
tile z a is locally asymptotically most powerful at asymptotic level a forHg against 

Uf >0 UcreR+ {Pe,cr,|;/i } • 

The advantage of the test statistic (3.14) compared to (3.10) is that, irrespective of 
the underlying density g\ , its behavior does not depend on the definition of the scale pa- 
rameter. The case of a Gaussian reference density (/i =4>x), however, is slightly different 
due to the particular form of the score function <pf 1 ; see Section 3.3. 



3.2.2. Unspecified location 

We now turn to the case under which both f\ and the location are unspecified. Again, 
9 is to be replaced with some estimator, but additional care has to be taken about the 
asymptotic impact of this substitution. Still, from LeCam's third lemma, it follows that 
the impact, under P^l . gi , of an estimated 9 on A^ 3 (k) can be obtained from the 
asymptotic behavior of 

which is asymptotically normal with asymptotic covariance matrix 

7£(/i) <&(/i.Si) 

where ^(/ljfli) := o-~ x (J gi (fi,gi) - KZ gi (fi,gi)). Clearly, this covariance 5 gi {fi,gx) 
vanishes iff k — K gi (/i, <7i), which, for g\ = f\, coincides with k(/i). 

Assuming that an estimate K^ n '{fi) such that — K gi (fi,gi) = op(l) un- 

der Pg"^ . Sl exists, Aj"? 3 (k^(/i)) is asymptotically equivalent to Aj". 3 (/t(/i)) un- 
der Pgl x , and asymptotically uncorrelated with (9, a, 0) and A^ 2 f) 0) (hence, 
asymptotically insensitive (in probability) to root-n perturbations of both 9 and a) un- 
der Pg^ojfli- :t fon ° ws lrom Section 3.2.1 that K^°(h,6,a) defined in (3.12) is such 
an estimator. The same reasoning as in Section 3.2.1 implies that this still holds when 
substituting, in A^™? 3 (k), any estimators Q# and <r# satisfying (CI) and (C2) for 6 and a. 
Finally, Lemma A. 5 in the Appendix ensures that A^". 3 (ft'"' (/i,##,<7#),##,<7#) can be 
substituted for Ay™. 3 (re si (/i,<7i),0,cr). We thus have shown the following result. 



Proposition 3.3. Let fx € F{° and gi 6 J 7 ^ . Then: 
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(i) t)" )o (0 # ,ct # ) =T<£ )o (6,& # ) +o P (l) = T { f " )o (8,a) +o P (l) is asymptotically nor- 
mal, with mean zero under Pi"i n .„ , mean ( 3. 15) under P„ 1/9 and variance 
one under both. 

(ii) TTie sequence of tests rejecting the null hypothesis : = U 9l e.7^ UeeR^e"gi °/ 
symmetry (with unspecified location 8, unspecified scale a and unspecified stan- 
dardized density g±) whenever T^" (##,<7#) exceeds the (1 — a) standard normal 
quantile z a is locally asymptotically most powerful at asymptotic level a for 
against (J € >0 UeeR IUr+ {PfcUji 1 ■ 

This test is based on the same test statistic as the specified-location test of 

Proposition 3.2, except that the (here unspecified) location 8 is replaced by an estima- 
tor The local powers of the two tests coincide: asymptotically, again, there is no loss 
of efficiency due to the non-specification of 0. 

3.3. Pseudo-Gaussian tests 

Particularizing the reference density f\ as the standard normal one 4>i in the tests of Sec- 
tions 3.2.1 and 3.2.2 in principle yields pseudo- Gaussian tests, based on the test statistics 
^i™ W or ^i™ -Due to the particular form of the Gaussian score function, however, 
the Gaussian statistic can be given a much simpler form. Indeed, I gi {4>i,g\) =I(4>\) = a 
does not depend on g\, and needs not be estimated, while J gi (<f>i,gi) = J{4>\) = 3a^2(si), 
so that K ffl (/i,<7i) is consistently estimated by '(0) / 'a 2 . This, after elementary com- 
putation, yields the test statistic 

n 

rwt ((9) := Y,{Xi - o)«Xi - ef - 3m< n) (6)), (3.16) 

where := y^(0) := mf ) {0) - 6m^ l) (8)112^ ) (0) + 9(m 2 n) (6>)) 3 . For this test statistic 
T( n ^(8), the asymptotic shift (3.15) under Pg" < 7n -i/2 T . gi now takes the form 

t[SM4 (31 ) - 9f4 (31)] [fx 6 (31 ) - 6^2 (.91 )M4 (ffi ) + 9^2 (.9i )] ~ 1/2 • 

This shift does not depend on a anymore, and still reduces to r-^/6/a under P^ j1 _i/ 2t . { ^ i 
(the same value as for T^(0, a), which confirms that optimality under Gaussian densities 
has been preserved); nor does it depend on the scale. 

The tests based on the asymptotically standard normal null distribution of T^™^ are 
optimal under Gaussian assumptions, but remain valid when those assumptions are vi- 
olated. Again, a simple Slutsky argument allows for replacing 8 (if unspecified) with 
any consistent estimator 8 without going through discretization; moreover, (3.16) does 
not depend on a. The tests based on T^{8) and T^(X^) both are closely related 
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to the traditional test of symmetry based on More precisely, under any ~Pg„ . gi , 
gi G {J 7 ^ ~)J'(j>i (note that the assumption g± £ J-^ = J- t j >1 implies that g\ has finite 
moments of order six), 

T (»)t(fl) = r Wt(lW) + 0p (i) = + 0p ( 1 ) j (3.i 7 ) 

where S^ 1 is the empirically standardized form (1.3) of bi (see (1-3)). 
Summing up, we thus have the following result. 

Proposition 3.4. Let g\ 6 Tfa and 8 = 6' + Op(n~ 1 / 2 ) ; recall that Hk{g\) '■= z k g\(z) dz 
stands for g\ 's moment of order k. Then: 

(i) rM+(<?) = T(")t(6») + o P (l) is asymptotically normal, with mean zero under 
F e,afi-,gi> mean 

r[5/i 4 (gi) - 9l4(gi)]/\Me(gi) - 6// 2 (ffi Wsi) +9^(.9i)] 1/2 

under P„ . „ and variance one under both. 

(ii) The sequence of tests rejecting the null hypothesis of symmetry (with specified 
location 9 ) T-tg := Ugie.F ^0™gi whenever T("^ (6) exceeds the (1 — a) standard 
normal quantile z a is locally asymptotically most powerful at asymptotic level a 
against U ?>0 U^r+^S,^}- 

(iii) The sequence of tests rejecting the null hypothesis of symmetry (with unspeci- 
fied location) := U ffl e^ UeeR^e"gi whenever T^ n ^{6) exceeds the (1 — a) 
standard normal quantile z a is locally asymptotically most powerful at asymptotic 
level a against U € >oUe GK U ff eR+{ P S,€tfi}- 

For the sake of completeness, we also provide (with the same notation) the following 
result on the asymptotic behavior of the (suboptimal) test based on m^™' (0) . Details are 
left to the reader. 

Proposition 3.5. Let g 1 6J^. Then, s[ n) := n 1 / 2 m| n) (0)/(r4™ ) (0)) 1/2 is asymptoti- 
cally normal, with mean zero under P# Jo-gi> m ean r[5/i4(gi) — 3/c(<7i)/i2(<7i)]/ Me/ 2 (&i) 
under P„ . and variance one under both. 

Under Gaussian densities (gi = </>i), the asymptotic shifts of T^°(9) (Proposi- 
tion 3.3(i)) and s[ n ^ (Proposition 3.5) are 16t/\/6 and 16r/\/l5, respectively; the asymp- 
totic relative efficiency of T^ n ^{9) with respect to is thus as high as 2.5 in the vicinity 
of Gaussian densities. This, which is not a small difference, confirms the suboptimality 
of (0)-based tests. 
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3.4. Laplace tests 

Replacing the Gaussian reference density (f>\ with the double-exponential one fc, we 
similarly obtain the Laplace tests. The assumption that f\ G J-° unfortunately rules 
out fc, since <pf c (z) = sign(z)/d is not differentiable, so that the construction of k^°(/i) 
in (3.12) does not apply for f\ = fc- Now, a direct construction is possible: l gi (fc,gi) 
indeed reduces to 2gi(0)/d - which is consistently estimated by I^°(fc) '■= 2<h(0)/d 
(where g\, e.g., is some kernel estimator of g\). Similarly, J gi {fc,9i) reduces to 
(2/d) \z\gi (z) dz - which is consistently estimated by J {n ^° (/ £ ) := (2/nd) £™=i \ z % 
<r#)|; the scaling constant d is easily computed, yielding d= l/(log2) w 1.44. Then, 

1 " 

^{fc)^J {n) %fc)/I (n) °{fc)^-^Y.\ Z ^^#)\ 

n 9iv J ) i=1 

is such that n^°{fc) — K (fc) = op(l) under .f c , as required. 

The Laplace tests are based on (specified 9) or T^°{9) (unspecified 9), where 

n 

if )o (0) := - Vsign(Z t (0,a # ))((Z t (g,a # )) 2 - ^'(/i)) (3.18) 

y/nrfWtf/:) t=1 

with 

(«) (n) n / 1 n \ 2 

These tests share with the Gaussian Fechner test (see [3]) the use of the score function 
z i — y sign(z)z 2 . The orthogonalization, however, differs, since the Fechner and Edgeworth 
families the tests were built on are different. The following proposition summarizes their 
properties; details are left to the reader. 

Proposition 3.6. Let g\ <E J~f c , 9 = 9 + Op(n -1 / 2 ) and denote by H\k\{gi) '■= 
J^° oa \z\ k gi(z)dz the absolute moment of order k of g%. Then, 



(i) T^ (8) = (9) + Op(l) is asymptotically normal, with mean zero under 

P 9,a,0 m > mean 

T[4» w (g 1 )-2tf l ,( gi )/g 1 (p)] 



[ M4 (3i)-2 M2 (. 9l ) / i| 1 |( 5l )/ 3l (0) + Ai 2 1| ( 5l )/(« ? i(0))2]i/2 



under P„ ,,, and variance one under both. 



,<r,n 1 / 2 r;^i 



(ii) The sequence of tests rejecting the null hypothesis of symmetry (with specified 
location 9) H^ 1 ' := U ffl e.?> ^H^gi whenever T^ (9) exceeds the (1 — a) standard 
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normal quantile z a is locally asymptotically most powerful at asymptotic level a 

against U c>0 IUr+ { P £,f ; / £ } ■ 
(iii) The sequence of tests rejecting the null hypothesis of symmetry (with unspeci- 
fied location) := U 9l e.F UeeR^e'gi whenever (9) exceeds the (1 — a) 
standard normal quantile z a is locally asymptotically most powerful at asymptotic 
level a against [j i>0 UemUam+i^l^U }• 

Comparing the asymptotic shifts of the pseudo- Gaussian tests and the Laplace ones 
yields asymptotic relative efficiency values; the asymptotic efficiency of tests based 
on T^"^ with respect to those based on is 1.76 in the vicinity of Gaussian densi- 

ties and 0.7 in the vicinity of double exponential ones. Finally, note that the empirical 
median X^,^ here provides a much more sensible estimator of 9 than the empirical 

mean XW; it has been used for 9 in the simulations of T^ l '°(§) in Section 4. 

4. Finite sample performances 

We performed a first simulation study on the basis of N = 5000 independent samples of 
size n = 100 from (2.2), with normal and double-exponential densities f\ and skewness 
parameter values £ = 0.1 and £ = 0.2. Each of those samples was subjected, at asymptotic 
level a = 5%, to the classical specified location test of skewness based on m^ 1 ' } (9) (i.e., 
on (1.1)), the (optimal) pseudo-Gaussian tests based on £>j n ' (i.e., on (1.3)) and the 
corresponding Laplace and logistic tests. For the sake of completeness, the two triples 
tests proposed by Randies et al. [13], which arc based on the signs of Xj + X, — 2Xk, 
l<i<j<k<n, are also included in this simulation study. Those tests, which are 
location-invariant, do not follow from any argument of group invariance and are not 
distribution- free. 

Rejection frequencies are reported in Table 1. 

Note that all tests considered here, except for Randies', are extremely conservative, and 
in most cases hardly reach the nominal 5% rejection frequency under the null. Randies' 
tests, on the other hand, significantly over-reject, which does not facilitate comparisons. 
Despite that, the tests based on b[ n) , T { c n)c ' {X { ™\) and T^J°{X^) exhibit excellent pe- 

formances, and largely outperform those based on (6) (although the latter requires 9 
to be known). 

The Edgeworth families considered throughout this paper, however, served as a theo- 
retical guideline in the construction of our Edgeworth testing procedures, and never were 
meant as an actual data generating process. One could argue that analyzing performances 
under alternatives of the Edgeworth type creates an unfair bias in favor of our methods. 
Therefore, we also generated N = 5000 independent samples of size n = 100 from the 
skew-normal SJ\f(X) and skew-t St(v, A) densities (with v = 2, v = 4 and v = 8 degrees 
of freedom) defined by Azzalini and Capitanio [1], for various values of their skewness 
coefficient A (A = implying symmetry); since the sign of A is not directly related to 



Optimal tests for symmetry 



1083 



Table 1. Rejection frequencies (out of N = 5000 replications), under various symmetric and 
skewed normal and double-exponential distributions from the Edgeworth families (2.2), with 
£ = 0, 0.1, 0.2, of the classical tests of skewness, based on m^"' (8) and 6^ , the Gaussian, Laplace 
and logistic tests, and the triples tests and of [13] 



Test 








SC(0 






t 















0.1 


0.2 





0.1 


0.2 




0.0372 


0.1136 


0.0996 


0.0306 


0.6938 


0.8722 


T {n) \6) 


0.0434 


0.7276 


0.9958 


0.0252 


0.4596 


0.6774 




0.0416 


0.6986 


0.9746 


0.0444 


0.7458 


0.8930 


r<"'°(«) 


0.0520 


0.5424 


0.9474 


0.0406 


0.9090 


0.9998 




0.0280 


0.4440 


0.8360 


0.0284 


0.8838 


0.9960 




0.0492 


0.7336 


0.9954 


0.0378 


0.8516 


0.9894 




0.0362 


0.6626 


0.9716 


0.0384 


0.8516 


0.9880 


rp(n) 
± R1 


0.0518 


0.6786 


0.9606 


0.0576 


0.9276 


0.9986 


rp(n) 
± R2 


0.0608 


0.6992 


0.9640 


0.0650 


0.9350 


0.9988 



that of £, we only performed two-sided tests. That class of skewed densities was chosen 
in view of its increasing popularity among practitioners. 

None of the tests considered in this simulation example are optimal in this Azzalini 
and Capitanio context. Inspection of Table 2 nevertheless reveals that the classical tests 
of skewness based on m^ l \9) and b± collapse under ti and £4, which have infinite 
sixth-order moments, and under the related St(2,X) and Si(4,A) densities. The same 
tests fail to achieve the 5% nominal level under the Student distribution with 8 degrees 
of freedom (despite finite sixth-order moments) and show weak performance under the 
St(8, A) density. Remark that the suboptimality of the test based on 171^(6), which, as 
a consequence of Proposition 3.1, may be considered as an artificial consequence of the 
choice of skewed families of the Edgeworth type, nevertheless also very neatly appears 
here. The triples tests behave uniformly well; note, however, their tendency to over- 
rejection, in particular under Student densities. 

5. Conclusions and perspectives 

We have derived the optimal tests for testing the hypothesis of symmetry within fam- 
ilies of skewed densities mimicking, in the Gaussian case, the type of local asymmetry 
observed in a central limit behaviour. The resulting tests were obtained for specified 
and unspecified locations, under specified or unspecified densities (satisfying appropriate 
assumptions). Local powers and asymptotic relative efficiencies are computed and finite- 
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Table 2. Rejection frequencies (out of N = 5000 replications), under various symmetric and related skew-normal and skew-t 
distributions [1] SJV(X) and St(u, A) (V = 2,4,8 and various A) of the classical tests of skewness, based on m ^(d) and ^ n) , the 
Gaussian, Laplace and logistic tests, and the triples tests T^™' and of [13] 

Test A 

SAf(X) St(2,X) St(4,X) St(8,X) 

0123024602460246 
0.0476 0.0482 0.0952 0.1936 0.0072 0.0106 0.0126 0.0144 0.0192 0.0190 0.0298 0.0436 0.0316 0.0608 0.1302 0.1754 
0.0374 0.0634 0.2988 0.5942 0.0046 0.0118 0.0182 0.0268 0.0144 0.0184 0.0586 0.1134 0.0260 0.1422 0.4078 0.5428 
0.0418 0.0616 0.3066 0.6130 0.0172 0.0232 0.0308 0.0396 0.0252 0.0302 0.0754 0.1298 0.0322 0.1604 0.4186 0.5484 t> 

0.0460 0.0690 0.2682 0.5406 0.0180 0.0414 0.0740 0.0988 0.0332 0.0406 0.1304 0.2514 0.0456 0.2054 0.5640 0.6906 O 

a 

0.0334 0.0472 0.2022 0.4736 0.0168 0.0288 0.0520 0.0678 0.0206 0.0302 0.0950 0.1798 0.0304 0.1518 0.4834 0.6260 | 
0.0468 0.0742 0.3542 0.7010 0.0154 0.0286 0.0496 0.0666 0.0236 0.0318 0.1082 0.2086 0.0342 0.2104 0.5846 0.7508 i 1 
0.0354 0.0568 0.2988 0.6426 0.0144 0.0256 0.0408 0.0492 0.0228 0.0276 0.0882 0.1694 0.0288 0.1712 0.5046 0.6696 ^ 
0.0540 0.0778 0.3602 0.7082 0.0618 0.1032 0.1798 0.2312 0.0508 0.0636 0.1842 0.3444 0.0530 0.2592 0.6766 0.8336 g? 

0.0598 0.0886 0.3812 0.7258 0.0656 0.1098 0.1882 0.2424 0.0556 0.0688 0.1938 0.3582 0.0598 0.2740 0.6940 0.8422 IF 

a 
s 
a. 

S3 

a 
^ . 
s 
a- 
a 
~ 
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sample performances are investigated in the context of classical skew-normal and skew-i 
distributions. 

Establishing the optimality properties of the traditional test based on the Pearson- 

(n) 

Fisher coefficient b\ was one of the objectives of this work, and wc show that this test is 
indeed optimal in the vicinity of Gaussian symmetry. Interestingly, its optimality holds 
for the specified-location as well as for the unspecified-location hypothesis of symmetry 
and is preserved if centering, in the computation of the test statistic, is based on robust 

- (n) 

root-n-consistent estimators of location rather than on the sample mean X 

Table 3 provides a summary of the various tests described throughout Section 3, along 
with their validity and optimality properties. 

These tests naturally extend into nonparamctric rank-based ones. The hypothesis of 
symmetry indeed enjoys strong group invariance features. The null hypothesis T-L^ of 
symmetry with respect to 9 is generated by the group Qi n \o of all transformations g h 
of K" such that sh(xi, x n ) := (h(x\ ),..., h(x n )), where \\m x ^± oa h(x) — ±oo, and 
x I— > h(x) is continuous, monotone increasing and skew-symmetric with respect to 9 (that 
is, satisfies h{9— z) — 9 = —(h(9 + z) — 9)). A maximal invariant for that group is known to 

be the vector (si(0), . . . , s„(0)), along with the vector (Rv\(6), . . . , R^^O)), where si{9) 

is the sign of Xi — 9 and R™\(9) the rank of \Xi — 9\ among \X± —6\,..., \X n — 0\. General 
results on scmi-paramctric efficiency [7] indicate that, in such a context, the expectation 
of the central sequence (t?) conditional on those signed ranks yields a version of the 
semi-parametrically efficient (at f\ and i9) central sequence. 

That rank-based approach is adopted in a companion paper [4]. For instance, the rank- 
based counterpart of the specified-^ test statistic of Proposition 3.6(ii) is the (strictly 
distribution-free, irrespective of any moment assumptions) van dcr Wacrdcn test based on 

rS(S) := ^m^" m ( )((* ( ^ )) - 3 )- 

where 7 (n) (0i) := E"=i ^Hf^^X^Hl^y)) 2 - 3) 2 and $ stands for the stan- 
dard normal distribution function. The unspecified-^ case under such approach, however, 
is considerably more delicate. 

Appendix 

A.l. Proof of Proposition 2.1 

The proof relies on [16], Lemma 1, which involves a set of six jointly sufficient conditions. 
Most of them readily follow from the form of local likelihoods and are left to the reader. 

1/2 

The most delicate one is the quadratic mean differentiability of (9, a, £) g s ' a , (x) , 
which we establish in the following lemma, where ge^^-j^x) is the density defined 
in (2.2). 



o 

oo 
— 



Table 3. A summary of the various test statistics considered throughout the paper, with reference to their definitions, the testing 
problem (specified or unspecified location - scale throughout remains unspecified) addressed, the standardized densities gi for 
which they are valid, the densities at which they are optimal and the proposition (lemma) where their asymptotic properties are 
described 



Test statistic 


Reference 


Location 9 


Validity (comments) 


Optimality (comments) 


Asymptotic properties 


o 1 [tf ) or m 3 [u ) 


(1.1) 


Specified 


Any g± with finite 


Suboptimal at the Gaussian 


Proposition 3.5 


(classical) 






6tli-order moments 






o 2 or t> 1 


(1.3) 


Unspecified 


Any gi with finite 


At the Gaussian, uniformly in 0,(j 


Proposition 3.4 


/ i • i ^ 

(classical) 






6th-order moments 


(=T K "(X ) + o P (l)) 


(3.17) 




(3.1) 


specified 


bpecmed /i 6 J~i 


At / CT , uniformly m it 


Proposition 3.1 


T ( f :\0#,a # ) 


(3.1) 


Unspecified 


Specified fi G T\ 


At fcr, uniformly in 9, a 


Proposition 3.1 


f^(9,a # ) 


(3.10) 


Specified 


Any g L G F* fl (3.9) 


At f a , uniformly in a (bad behavior 


Lemma 3.1 








under local alternatives and g\ 7^ /1) 






(3.14) 


Specified 


Any Sl e J? (3.13) 


At / CT , uniformly in cr (same as t|™ , 


Proposition 3.2 








except for its bad behavior under 
local alternatives and f/i 7^ /1 ) 






(3.14) 


Unspecified 


Any 31 <E J? (3.13) 


At f a , uniformly in 9, a 


Proposition 3.3 




(3.16) 


Specified 


Any 31 £ (implies 


At the Gaussian, uniformly in cr 


Proposition 3.4 


(pseudo-Gaussian) 






finite 6-order moments) 






T (n)t ((9) 


(3.16) 


Unspecified 


Any gi G Tfa (implies 


At the Gaussian, uniformly in 9, a 


Proposition 3.5 


(pseudo-Gaussian) 






finite 6-order moments) 






T ( c n)o {6) 


(3.18) 


Specified 


Any gi 6 Tj L (implies 


At the double exponential, 


Proposition 3.6 


(Laplace) 






finite 4-order moments) 


uniformly in a 




)o (8) 


(3.18) 


Unspecified 


Any gi 6 Tf c (implies 


At the double exponential, 


Proposition 3.6 


(Laplace) 






finite 4-order moments) 


uniformly in 9, a 





^ . 

s 

(Is 
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a I a 



-K(h))I[\x-6\<a\z* 



- sign^a- 1 /! f — - A > sign(-Oa|^|] 

-J[a;-fl<8ign(0^*l]}. 

1 -3/2 f l/2{ x-6\ fx -6" 
2 a h \ a 



—a 
2 



- 1 



.f*-/ [ "':"(T)*.(T)((v)'-«' 



Then, as r, s and t — > 0, 
(i) / {9g+t,<r+,,T i f 1 ( x )-9g+t,<r+,,o i f 1 ( x )-^ 

I / n 1/2 



dx = o 



W£, CT+s ,i ; / 1 Wle=o-^i 2 ,c ; / 1 (^)le=o)} 2 da ; = o(i). 



(iv) / <! ffo5 1<r +. lP; / 1 (af) ~ 9eJ, ;fM ~ I - s 



1/2 



D "9 l eXo- h { x ) 
\ D £9l / *, e ;f 1 ( x )k=o, 



dx 



Proof, (i) Decompose j {gl+ t ^ +s . r , h (x) - .9e+l ff+sA /i ( x ) ~ rD t9e+t,<r+s,t;fi( x )} 2 dx 
into ai + 2a2, where 



1/2 



1/2 



ai = 



. { (^/iN) 1/2 [i + ^ (")(« 2 - *(/i))] 1/2 - (^/iN) 



1/2 
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2 fl'\u) J 



and 

02= [ {((a + s )-Vi(«)) 1/2 -^(^ + S )- 1 / 2 4#7(« 2 - K (/i))j 2 (^ + s ) d ^ 
J«>|z*| I z // (w) J 

Since, for |x| < 1, (1 + x) 1 ' 2 = 1 + f (1 + \x)- 1 ' 2 for some A £ (0, 1), one easily obtains 
that 

= l(-+ s r 1/2 -^#r(- 2 --(/i)) 



4 



|u|<|z* 



A 1/a («) 



x((l + Ar^ A («)(u 2 -K(/i))) 1/2 -1)S> (a + s)du. 



For |u| < 1, one has (1 — (1 + Alt) 1 I 2 ) 2 < 2j3j, and the integrand is bounded by 

2 i r(" - K (/i)J 



which is squarc-integrable; Lebcsguc's dominated convergence theorem thus implies 
that a\ is o(r 2 ). Turning to a2, we have that a2 < C((cr + s)~ 1 a2i + 022), where 

f ^* 2 f ( f (^) \ 2 
a 2 i := / fi(u)du and a 22 := — / ( — T75 ) (w 2 - k(/i)) 2 du. 

J«>iz*i 4 y„>|«.| V/i /2 («)/ 

The definition of T\ implies that 021 = 0((z*)~P), hence that a 2 i = o(r 2 ) if r(z*) /3 / 2 — > 00 
as r — > 0. This latter condition holds, since 4>f 1 (z) = o(z' 9 / 2 ~ 2 ) and since the definition 
of z* entails that 

1 -wyyva ( Z * 2 - K (A)) 

^ ' z *(p/2-2) z *2 

An application of Lebcsguc's dominated convergence theorem again yields 022 = o(r 2 ). 

(ii) This is a particular case of Lemma A.l in [6] (here in a simpler univariate context). 

(iii) The fact that D^g^^ ^ .,(a;)|^ = o is square-integrable implies that 

\\ D i9e+t,a,uS x ^ - D t9e / ,*,e;f 1 ( x )\t=oh*=o(l) 

as t tends to zero. Define fi ;e xp(x) := fi{e x ) and (fl{^ p (x))' := |/i 2 (e x )/i(e x )e x . For 
the perturbation of a, we have 



J {( B e5^a+« > c;A( a )le^- D e5»5,e ; /i( :B )le=o)} 2da! 
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-3/2 



-2a 



o 



~) 3 2 (/ 1 1 ; / exp)'(lnW-ln(l + ^)) - (/iYcx P )'( ln W) 



((i + J) 1 2 (/lYexp)' (M*) - ln (l + J) ) + (O'OK*)) 



2 

dz 



<C( Cl +c 2 ), 
where 



Cl 

and 

C2 = 



" (e<»™-M l+ ./->) (/*£,)' (« - ln(l + i) ) - e^-C/^Hu)) ' d« 

_ (e-(V2)(-Mi +S / CT )) tf/^y ( u _ + £^ _ e-C^^)'(«)) ' d«. 

Now, both e _ ( 1 / 2 ) ll (/ 1 1 ; {x P )'(w) and e^ 3 / 2 ^/^ ar <3 square-integrable since fi&Fi- 
Therefore, quadratic mean continuity implies that c\ and C2 are o(l) as s — > 0. 
(iv) The left-hand side in (iv) is bounded by C(&i + 62 + ^3), where 



2 



and 



&3 = y M J3 eff0+t,«r+»,e ; / 1 ( a O - D ^9 1 e / ^. Jl {x))} 2 dx. 
The result then follows from (i)— (hi). □ 

A. 2. Asymptotic linearity 

A. 2.1. Asymptotic linearity 0/ A^™' 3 

The asymptotic linearity of A^™? 3 is required in the construction of the optimal parametric 
test of Section 3.2.2. Note that the proof below needs uniform local asymptotic normality 
in 9 and a only. 

Proposition A. 1. Let fx € T\ and gi G J 7 f 1 - Then, under . gi , as ?i— >oo ; 
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(i) A^3(e + n-V%a,0) = Ag 3 (^a ) 0)-t ( 7- 1 (J gi (/i,5i)-«(/ 1 )X gi (/i,ffi)) + op(l) 
for allteR. 

(ii) A ( £] 3 (e,o- + n-V 2 s,0) = ^ { f %(0, er, 0) + o P (l) for all s G R. 
Proof. Define 

n 

<»3 (*, *) := n- 1 /^ ^ ^ (Z W (fl> CT)) (z (") (0 , a)) \ 

i=l 

Letting A^u) := ^ (G*r+( u )) ( G i+( u )) 2 with G*i+(z) := (2Gi(z) - > 0], note that 

n 

Df^{6^Q)=n- 1 '^ s M K fAGx + {\z\ n \e,a)\)). 



Writing 0W for tf + n- 1 /^, z| n) for } (6>,a), sf } for sign(zl™ ) ), for zj n) (6> (n) ,CT) 



(n) 



and s-!^ for sign(Z^), let us show that 



" 1/2 E4 ) ^ 1 ( G i+(i^£ ) i))-"" 1/2 E s l" ) ^( G i+(i^" ) i))+ to " 1 ^(/i^i) 



i=l 



is op(l); the proofs of (ii) and that of 

n n 

n- 1/a E^(^5) = »" 1/a E^i(^ B V^-%(/i,ffi)+op(l) 

i=l i=l 

follow along the same lines and are therefore left to the reader. Let 



+ jr /l (i- T ]i( (i_ 7 ]_ tt ]j 



1 2 

7 <u£ 7 



i 



2 2 



2 1 

1 -y < ^ 1 -y 



Continuity of u H> Kf^u) implies continuity of u i-> K^(u) on the interval ]0, 1[. More- 
over, since this function is compactly supported, it is bounded, for any (sufficiently large) 
/ € No, by the monotone increasing function u \— > (it). Let Eo denote expectation under 
P^l Q . One easily shows that D^Pl 3 (9,a) decomposes into £)^"' m ' + £)^ n,m ^ — R^ 1 '™' + 

R (n,rn) +R (rn) jWfagJ .— K^(u)(j) gi (G^ 1 (u)) du, 

n n 
i = l i=l 

-n^E a [s i;n K^(G 1+ (\Z^\m 
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4».0 =n -V2^ s . (ir/i(Gl+( | Z (n) |)) _ ir 0) (Gi+(|z (n) |)))) 
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i 4».0 =n -i/ a ^^ n( ^ (Gl+( |^)| )) _^(0 (Gl+( | Z W| ))) 



and 



R$ ) =to-- 1 (J gi (fu9i)-J<; l Mh,9i)) 



In order to conclude, we prove that D^' and D^' are op(l) under Pg ^ , as n — > oo, 

for fixed Z, and that Ri' 1 ' , R^ 1 ' 1 ^ and -R3 are op(l) under the same sequence of hypothe- 
ses, as I — > 00, uniformly in n. For the sake of convenience, those three results are treated 
separately (Lemmas A.2-A.4). 
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Lemma A. 2. For any fixed I, D^ 1 ' 1 ' = Op(l) as n —> 00, under Pgl 
Lemma A. 3. For any fixed I, D^ 1 ' 1 ^ = op(l) as n—> 00, under Pi „ n .„ . 

Lemma A. 4. (i) Under Pg"^ . fll , iJj™ 1 =op(l) as Z— >oo, uniformly in n. 

(ii) i?2™' i ' ) = op(l) as / — > 00, under Pg"J ^/or n sufficiently large), uniformly in n. 

(iii) is o(l) asZ—> 00. 



Proof of Lemma A. 2. Consider the i.i.d. variables 



T- 



(n,l) 



Si ., n K^(G 1+ (\Z^\))- Si K^(G 1+ (\Zr\)). 



An) 



(I) I 



An), 



One easily verifies that D[ n ' l) = n" 1 / 2 £" =1 (i; (rM) - E [i; ( "' ]). Writing Var for vari- 



ances under P^ - 31 , we have that 



g (2 f».0_ Eo ^.O ]) 



< n 1 Var 



(n,Z)x2n 



and it only remains to show that 
Vol(Tj n -y}=E {(s l;n ^ 
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as n — > oo. Now, 

( Si . >n K^(G 1+ (\Z^\))- Si K^(G 1+ (\Z^\))) 2 

= ( Si ., n K®(G 1+ (\Z$\)) - atsn K®{G 1+ (\Zy>\)) + Si ., n K^(G 1+ (\zl n) \)) 

- Si K^(G 1+ (\zl n) \))) 2 

< 2(K${G 1+ (\Z$\)) ~ Kf{G 1+ {\Z^\))f + 2(K^(G 1+ (\Z^\))f(s i;n - Sl ) 2 . 

Because u^K$(u) is continuous and \Z$ - z\ n) \ is o P (l), K%(G 1+ (\Z$\)) - 

Kfr \Gi+(\Z\ ; I)) also is op(l). Moreover, since is bounded, this convergence to 

zero also holds in quadratic mean. Similarly, K^{Gi + {\Z[ n '\)){s^ n — s.;) = op(l) since 

K®(Gi+(\Zi\)) is bounded and |sj ;n — Sj| is op(l). Finally, both s i;n and Sj are 
bounded, implying that this convergence to zero also holds in quadratic mean. □ 

Proof of Lemma A.3. Let B[ n ' l) := n~ 1/2 £™ =1 s l Kf i {Gi+{\Z i i n) \)). Asre^oo, under 
P^l . gi , Bi' is asymptotically Af(0, E[(AT $ (C/)) 2 ]), where [/ stands for a random vari- 
able uniformly distributed over the unit interval. Also, letting B% :— n" 1 / 2 s «:« x 
K^(G 1+ \Z^\), it follows from ULAN that B { 2 n ' l) - tcj- 1 j£\fi,gi) is asymptotically 
Af(0,E[(irJ?V)) 2 ]) as n -> oo, under pg, 0;si • Since D^"'^ = B^ l) - B[ n ' l) -E [B^ l) ] = 
o P (l), we have that B { 2 n ' l) -E [B { 2 n ' l) ] is asymptotically Af(0, E[(K { f [ } (U)) 2 }) as n->oo, 
under P^,o ;9l - Therefore, still asu->co, D^™' = E [S^ n ' ] -ta^J^ifugi) = o(l). □ 

Proof of Lemma A. 4. (i) We have that 
Eop^) 2 ] < CEoKi^^ 

r OO 

= C/ (#*(«) -l$(t0) 2 d«; 



for any uG]0,l[, (u) converges to Kf 1 {u) and the integrand is bounded (uniformly 
in I) by 4(AT/ 1 (u)) 2 , which is integrable on ]0, 1[. Lebesgue's dominated convergence 
theorem thus implies that Eo[(-R^"''') 2 ] = o(l) as I — > 00, uniformly in n. 

(ii) The claim here is the same as in (i), with zj^ replacing zf 1 ' . Accordingly, (ii) holds 
under P;?;-, „ . That it also holds under Pi"2n „ follows from Lemma 3.5 in [91. 



(hi) Note that 



\J g Ah,9i)-j£(fi,9i)\ 



^ gi (G^(u))(K^(u) - K fl (u))du 



<l( 9l ) I ((K^(u)-K fl (u))f 



dit, 
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where the integrand is bounded by 4(K ^(u)) 2 , which is square-integrable. Pointwise 

convergence of K^(u) to Kf 1 (u) implies that J^g 1 {fi,gi) — Jg} (/i)<7i) =o(l) as I — > oo. 
The result follows. □ 

A. 2.2. Substitution of A%] 3 ( K ( n > (fx, # ,<7 # ),0 # , <r#) for ^(K gi {h,g x ),6,o-). 

(n) 

Lemma A. 5. Let fi G J 7 ! and gi £ J 7 ^ . Then, under P e a . gi : 

(i) Ag 3 ( Ksi (/ 1 , 5l ),%^ # )-Ag 3 (K sl (/ 1 , 5l )^^) = op(l). 

(ii) AW( K W°(/J #1 %)j # ,%)-AW(/ t9l (/ 1)5l ) ) fl #) %)=op(l). 
(hi) Ag3( K (") o (/ 1 ,%,a # ),%,a # )-Ag 3 ( K3l (/ 1 ,.g 1 ),0,a)=op(l). 

Proof. Part (i) is a direct consequence of Proposition A.l. The left-hand side in (ii) can 
be written as 



2f> x If' := (k ( " )o (/i,%,(7 # ) - x n- 1/2 X)^(^i H**.**)) (A- 1 ) 

i=i 

where T^" 5 is o P (l). Now, ULAN implies that 

TW=n- 1 /^ 0/i (zW(^ ( x)) + ( ( x-% l( / 1 , 5l ),O)n^ - +0p( l), 

(A.2) 

(n) 

as n — > oo under P g a . Hence, the central limit theorem and the root-n-consistency 

of 0# and entail that (A.2) is Op(l); the result follows. As for (hi), it is a direct 
consequence of (i) and (ii). □ 
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